Entropy Bound for the TM Electromagnetic Field in the Half 

Einstein Universe 



I. Brevik*, R. Herikstad and S. Skriudalen 

Department of Energy and Process Engineering, 
Norwegian University of Science and Technology, N-7491 Trondheim, Norway 

(Dated: February 1, 2008) 

Abstract 

An explicit calculation is given of the entropy /energy ratio for the TM modes of the electromag- 
netic field in the half Einstein universe. This geometry provides a mathematically convenient and 
physically instructive example of how the electromagnetic and thermodynamic quantities behave 
as a function of the nondimensional parameter 5 = l/27raT, a being the scale factor and T the 
temperature. On physical grounds (related to the relaxation time), it is the case of small (5's that is 
pertinent to thermodynamics. We find that as long as 6 is small, the entropy/energy ratio behaves 
in the same way as for the TE modes. The entropy is thus bounded. The present kind of formalism 
makes it convenient to study also the influence from frequency dispersion. We discuss an example 
where a sharp cutoff dispersion relation can in principle truncate the electromagnetic oscillations 
in the Einstein cavity such that only the lowest mode survives. 
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I. INTRODUCTION 



Valuable insight can often be obtained by analysing models that are mathematically 
simple but nevertheless able to show the essentials of the physical properties of a system. 
The purpose of the present paper is to present an explicit calculation of the energy E'^^ 
associated with the transverse modes (TM) of the electromagnetic field in a spherically 
symmetric volume endowed with the metric of the half Einstein universe, and to combine 
this with the analogous transverse electric mode result E'^^ calculated earlier in order 
to obtain the total electromagnetic energy E = E'^^^ + E'^^ in the volume. 

The calculations will be given for a finite temperature T. One of our motivations is to 

n 

show, in the case of the TM modes as well as in the previous case of TE modes that the 
ratio of entropy 5*^*^ to thermal energy E'^'^^ is limited in the following sense: 

at high temperatures. Here 3 = 1/T. (Geometric units are used.) The half Einstein uni- 

nnnri 

verse, having received considerable attention in the past (see, for instance, Refs. [11,12,131,141), 
is both mathematically convenient to handle and provides an instructive physical example 
of how the entropy content is limited. The situation to be considered here thus falls into line 



with the current discussion in general about how entropy anc 



SB. 



energy behave both quantum 



In particular, a topic that 



mechanically and thermally in conformal field theories 
las attracted particular interest is the Verlinde bound for the ratio of entropy to energy 
sl. To what extent the entropy bound is realized in physical situations is usually not clear 
beforehand; it has to be investigated by concrete calculations in each case. Previous calcula- 
;ions on the entropy bound have usually been done in the high temperature approximation 
6|, lZI]. We shall also be assuming that T is high. This is actually a constraint that follows 
from thermodynamics: in order for thermal equilibrium to be established, T has to be much 
higher than the inverse transit time for light across the linear dimension of the system. This 
point is discussed in more detail in section II. 

As mentioned, the case of transverse electric (TE) modes in the half Einstein universe 
was considered in Ref. and the bound of Eq. ([T]) was verified in that case. The analogous 
result in the TM case is thus not unexpected, although the combined electromagnetic and 
thermodynamic behavior is complicated and can hardly be understood merely by inspection, 
without calculation. As we shall see, at higher orders in 6 (see Eq. ([3]) below), the TM and 



TE modes behave somewhat surprisingly, in the sense that the terms involving and 
in the energy do not cancel each other. This contrasts the behaviour known from Casimir 
theory in flat space, where the TE and TM modes compensate each other with respect to the 
mentioned terms and only the term proportional to T survives in E to the leading order (cf., 
for instance, equation (8.39) in [8j or equation (4.44) in {q]). The reason for this behaviour 
has to be related to the geometry of the half Einstein universe. 

An earlier version of the present paper lOl actually gave an impetus to the recent inves- 
tigation of Dowker on general spacetimes . Dowker's analysis is much more general than 
the concrete example on Maxwell fields considered in the present paper. We think, however, 
that there are definite advantages of going into concrete detail as we are doing here: for one 
thing, it becomes easy to analyse the influence from dispersion in the wall material. Usually 
in field theory, one assumes that the walls have infinity conductivity for all frequencies. This 
assumption is unphysical, as there is always a frequency dispersion present. Of importance 
in a Casimir context, is the behaviour of the permittivity e{i() as function of positive imag- 
inary frequency (. For a dielectric, the susceptibility ^(^C) — 1 is in essence proportional to 
(1 + ('^/uJq)^^, where ujq is the (dominant) resonance frequency. Typically, cuq ~ 1.5 x 10^^ 
rad/s. Thus ujq acts as a soft high-energy cutoff. When the frequencies are much higher 
than the resonance frequency one has ^(^C) " 1 = ^p/C^ fo^' bodies, metals or dielectrics, 
where Up is the plasma frequency. Taking gold as an example, Up = 1.37 x 10^^ rad/s. For 
very high frequencies the photons do not "see" the metal, and the number of modes in the 
cavity is therefore truncated. We shall return to this point in the final section. 

Let us also make some remarks on conformal invariance. The electromagnetic field is 
known to be conformally invariant in D = 4 spacetime dimensions. Thus the trace of the 
electromagnetic energy-momentum tensor is zero when D=4. In higher dimensions the sit- 
uation becomes more complicated. The Casimir energy for D > 4 was first calculated by 



Ambj0rn and Wolfram 



12j. We refer also to the recent analysis of Alnes et al. on the 



electromagnetic field between two parallel hyperplanes in higher dimensions, considering 
both metallic and MIT boundary conditions [13]. Using the axial gauge, the pressure be- 
tween the plates was found to be constant, while the energy density was found to diverge 
at the boundaries. This peculiar behaviour is a direct consequence of the lack of conformal 
invariance when D > A. 

Another related point worth noticing is the growing interest in higher derivative elec- 
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trodynamics, such as the Lagrangian form {F + 1/F) studied by Novello et al. [l^ (here 
F = F^jjF^^'^). One of the motivations of this kind of approach, among other things, is 
to model a phase of cosmic current acceleration. It would be of interest to apply such a 
theory also to the half Einstein universe, although we will not go further with this topic 
here. The approach has strong similarities with the {R+ theory in gravity, which has 
)een thoroughly studied in the recent past (see, for instance, papers of Nojiri and Odintsov 
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II. THE HALF EINSTEIN UNIVERSE 

The full static Einstein universe is of general physical interest as it is conformally equiv- 
alent to all closed FRW metrics. The Einstein metric is 

ds^ = -dt^ + {dx^ + sin^ xd^"^), (2) 

where x ^ [0, tt], 6* G [0, tt], and G [0, 2n]. The scale factor a is related to the cosmological 
constant A through a = A^^^^ = (AnGpo)^^^'^, where po is the energy density of matter 
(dust). The vacuum part of the energy density is A/SttG; the pressure is p = —A/SnG. 

The /ia// Einstein universe is characterized by x being restricted to one half of the previous 
interval, x ^ [0,7r/2]. A general point worth noticing in connection with this universe is 
;hat it is conformally related to anti-de Sitter space and is of relevance to supersymmetry 



The universe can be envisioned as a three-dimensional spherical volume spanned by 
the 'radius' x ^^id the angular coordinates 6 and 0, closed by a two-dimensional spherical 
surface at x = ''i"/2- We will take this surface to be perfectly conducting, this being an 
analogy to the Dirichlet boundary in the case of a scalar field. 

As there are two dimensional parameters in the problem, viz. f3 and a, it becomes natural 
to introduce 



(3) 

2™ ^ ^ 

as a nondimensional parameter. At high temperatures 5 <^ 1, and can thus be used as 
expansion parameter in a perturbative analysis. As anticipated above, it is the case of small 
(5's that is of thermodynamic interest. This can be seen from the following physical argument: 
In order to apply thermodynamic formalism to a fluctuating quasi-classical system, the 
temperature has to be sufficiently high to satisfy the condition T ^ 1/r, r being the 



relaxation time 



16| . In our case we may take r to be of the same order as the transit time 



for hght across a distance a, i.e., t a. [This is essentially the same kind of situation as 
experienced when a narrow beam of light impinges upon a liquid surface and makes it bulge 
outward; cf. the classic light pressure experiment of Ashkin and Dziedzic and also the 
theoretical discussion of it in Ref. 18|. In that case, the relaxation time is of the same order 



as the transit time for sound to traverse the beam.] Thus, in the present case T ^ 1/a, 
which means that 

5 < 1 (4) 

is the condition for applicability of thermodynamics. 

It is of interest at this point to make a comparison with the Casimir effect. Let us assume 
the usual Casimir configuration, namely two parallel metal plates separated by a gap a. At 
finite temperature we are dealing with discrete Matsubara frequencies, Cm = 2TTmT, with 
m > an integer. It turns out that the most important contributions to the Casimir force 
come from the frequency region (mO- ~ 1, or m ~ l/(27raT) (cf. the discussion on this 



point in 



19(1). That is, m ~ 5. Then, we may interpret the parameter 6 physically as the 



magnitude of the most important Matsubara numbers that occur in the analogous Casimir 
effect. This correspondence appears to be physically reasonable, though not trivial. 



III. GOVERNING EQUATION. THE TM MODES 

We take an orthonormal basis, {adx, cl sin x^O, a sin x sin 6d(j)), and split off the time factor 
as e"*"^*. From Maxwell's equations in curvilinear space we can derive the governing equation 
for (or H^). Denoting collectively these field components by X, and writing the angular 
contributions as spherical harmonics, X{x,0,(f)) = X {x)Yim{0 , <f)) , we obtain pj 

— f sin^ xX) + (uja? sin^ xX - 1(1 + 1)X = 0. (5) 
dx ^ ' 

The solution to this equation is 

X(xsin'-ixCS'^(cosx), (6) 

where n is an integer and C^^P are the Gegenbauer polynomials [2^ satisfying the differ- 
ential equation 

(1 - a;2)C^-)"(x) - {2a + l)x Cp^''^' (x) + p{p + 2a)C^''\x) = (7) 
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for p > 1. Inserting X into Eq. ([5]) we obtain the eigenfrequencies 

n + 1 

uJn = n>l>l. (8) 

To avoid infinities at the origin x = 0, we must have n — / > jj]. 

So far, the boundary conditions have not been considered. The transverse magnetic 
modes are subject to the boundary condition 

5^(sinxH^)=0, x=\. (9) 

where is the magnetic field component transverse to the radius x- From Maxwell's 
equations in orthonormal basis ^ we find, when comparing with the general solutions of 
the governing equation ([5]), that the magnitude if^ of the vector must be of the form 

ifxocsin'xC^+/)(cosx). (10) 

Thus we have from Eq. ^ 

9^{sin'+^xC^i-/^(cosx)} = (11) 

at the boundary. Observing the recursion relation for the derivatives of the Gegenbauer 
polynomials we get 

9.CW(0) = (n + 2a-l)Ct\(0) (12) 
at the boundary. The condition ffTTl) now yields 

= 0. (13) 

As the Gegenbauer polynomials vanish for odd m, {n — /) in Eq. (fT3l) must be even. 
It follows that for n even (odd), / must be even (odd) and the degeneracies are 

g^^^= (2/ + l) = ^(n + 2), n even, (14) 

^i°)= (2/ + l) = ^(n + l), n odd. (15) 

1=1,3,5,... ^ 

This leads to the following expression for the logarithm of the partition function Z'^^^ for 
the TM modes: 



InZ™ = - X](2n - l)(n + 1) In [1 - e-^-^^^+i)-^] 

oo 

-^(2n-l)nln[l-e-^™'']. (16) 



n=l 



n=l 



From this we identify two types of sums: 

oo 

G^") = ^n'^ln[l-e-2-(2«+i)5], a = 0,1,2, (17) 

n=l 

and 

oo 

= ^(2n- l)"ln[l -e"^""^], a = 1,2, (18) 

n=l 

so that 

InZ™ = -2G(2' - + - i(if(^' + H^'^). (19) 
We expand the logarithm in Eq. (1171) and take the derivative with respect to S: 

n=l fc=l 

1 r °° 
= - ds(27r5)-'^ r(s)C(s) J2 ii"(2n + 1)'"', (20) 

C(s) being Riemann's zeta function. In the last step above we made use of the relation 

e-^ = -^ / dsx-^Tis), (21) 

ZTTZ Jc 

where the integration contour is a line parallel to the imaginary axis at a sufficiently large 
value of 3fts. 

As we need to work out the sum over n it will prove worthwhile to derive a general 
expression for the integral in Eq. fl20|) . To this end we change the summation variable, 



En"(2n + l)^-^= E k'-'- (22) 

n=l fc=3,5,... \ ^ / 

The terms in brackets can be expressed as a binomial series. We add and subtract the k = 1 
term as well as the even terms, and insert the result into Eq. ( !20|) to get 

(9G(°) _ 1 
05 ~ ¥'i Jc 



ds(27r5)-^r(s)C(s) 



1=0 







(:) 


(-!)'{ 







\l—s+a—l 



)C(s-l-a + /)-l}. (23) 



When a = this expression becomes 



d6 



\ ds(2vr5)-T(s)C(s) {(1 - 2i-)C(s - 1) - l} . (24) 



The two terms in the integrand have poles for s = 0, 2 and s = 0, 1 respectively. There are 
thus three poles in all. 
When a — 1 we have 

dS 2i Jc 

The poles are at s = 1, 3 in the first term, and at s = 0, 2 in the last term. 
The final integral is for a = 2, 



1 /^ds(27r5)-T(s)C(s) {(1 - 2^-)C(s - 2) - (1 - 2^-)C(s - 1)} . (25) 



85 



= l|^ds(27r5)-r(s)C(s) (26) 
X {(1 - 2'-%{s - 3) - 2(1 - 22-)C(s - 2) + (1 - 2'-%{s - 1)} , 



with poles at s = 0, 4, s = 1, 3 and s = 0, 2 in the three terms respectively. Calculating all 
residues and collecting terms we find 

05 ~ 2452 (5 + 12 ' ^ ' 

dG^^^ _ C(3) _ J_ n_ 

~d5~~8^ 48^2 + 245^24' ^ ^ 

9^(2) TT C(3) n 1 vr 



85 9605^ Sn^^ ' 96^2 245 160' ^^^^ 
We turn now to H^"\ following the same steps as for G^"'\ First, we express the derivative 
in the form of an integral, 

= -j^ d.(47r(5)-T(.)C(.) E n'-\2n - 1)". (30) 

Again using the binomial series in the n sum we arrive at the following generic expression: 

2 



f ds(47r(5)-T(s)C(s) 



a 

1=0 



(a 
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-l)'2"-'C(s- 1 -a + O- (31) 



For q; = 1 we here get 

= -J^ ds(47r5)- r(s)C(s) {2C(s - 2) - C(s - 1)} , (32) 

with poles at s = 1, 3 in the first term and s — 0, 1, 2 in the second term. Similarly, for 
a = 2 we get 

ai/(2) 2 
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= - / ds(47r5)- r(s)C(s) {4C(s - 3) - 4C(s - 2) + C(s - 1)} , (33) 
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with poles ats = 0,4, s = l,3 and s = 0, 1, 2 in the first, second and third term respectively. 
Calculating all residues we obtain 

dH^^'^ C(3) vr 1 71 



dS 2452 '35 6^ 



(34) 



06 24054 27v^^ 246^ 65 10' ^ ^ 

Inserting the various terms into Eq. f|T9l) we obtain the following expression for the energy 
^TM ^ -d/d(3\nZ^^' of the TM modes: 

~ 24054 47r253 245^ + 245 240 ' ^ ' 

We similarly calculate the free energy F'^'^^ = — InZ^^: 

-"7205^+8^ + 245 + 24^^^" W^' ^^^^ 
and finally the entropy ^™ = p^dF^^^df]: 

- 1805^" 8^" 125 "24^^^+ 24- ^^^^ 



To leading order in 5 this yields 

^TM 4 



/5, 5 < 1. (39) 



3' 

Thus, in the high temperature limit the entropy of the TM modes is bounded, just as for 
the TE modes; cf. Eq. ([T]). 



IV. COMPARISON WITH THE TE MODES. THE TOTAL FIELD QUANTITIES 

The equality of the entropy/energy ratios for the TM and TE modes at high temperatures 
- meaning physically, as we have seen - that the temperature T is much higher than the 
inverse relaxation time 1/r - is as we might expect in view of the separability of the TM 



and TE modes in spherical geometry (cf., for instance, Ref. [2l|]). But according to our 
calculations there are differences between these modes as regards higher order terms in 
5 = j3/27ia. The expressions pertaining to the TE modes are 



.TE ^ C(3) TT , 13 4l7r 



24054 47r253 2452 + 245 240' 
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^^TP TT C(3) VT 13 , ^ 4l7r ^ . . 

3C(3) vr 13 , ^ 13 
* =ip-s¥-m-24'"*+24- f^^' 

These results were obtained in [l| via the same method as above, and also, as an independent 
check, via use of the Euler-Maclaurin sum formula. Adding the contributions from the TM 
and TE modes we obtain the following total field quantities, when expressed in conventional 
units, 

E = ^a'T^ - 2C(3)a^T^ - —aT^ + -T - i^, (43) 
15 ' 6 12 240a' ^ ^ 

F = -—a'T^ + C(3)a2T3 + -aT^ - -Tln(27raT) - (44) 
45 ' 6 12 ^ ^ 240™' ^ ^ 

S = ^a^T^ - 3C(3)a2T2 " Y «^ + i ^^(^TraT) + 1^. (45) 

Surprisingly enough, the TE and TM contributions to the and terms in the expression 
( 143|1 for E do not cancel out. 



SUMMARY AND FURTHER DISCUSSION 

Let us first summarize a couple of points: 

• We have assumed the parameter 5 = /3/27ia = l/(27raT) to be small. This is in 
accordance with the requirement of classical thermodynamics: under equilibrium conditions 
T has to be much larger than the inverse relaxation time [l^. When comparing with the 
Casimir effect between two metal plates, S may be given a physical interpretation as the 
magnitude of the most dominant Matsubara numbers jiol . 

• The most striking result of the above calculation is that the TE and TM modes do not 
compensate each other to orders and in the expression ( H3i) for the total energy. One 
might expect beforehand that the mentioned compensation should take place here as well 
as in the known case of flat space, considered earlier in connection with the Casimir effect 

We attribute the non- compensation to the properties of the Einstein metric. The 
formalism is generally too complicated to be transparent beforehand (it may be noted here 
that the degeneracies of the TM modes as given in Eqs. ( |T^ and ( fT5|) are complementary 
to those holding for the TE modes [l|). 

• We shall consider below some aspects related to dispersion. As a preliminary step, let 
us give first a brief account of the essence of the formalism for calculating the total energy 
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associated with the individual TE and TM field oscillation modes, characterized by the 

n 

numbers n and / ^1. 

n 

Consider first the TE modes. As in ^] it is convenient to change the meaning of n, such 
that n runs from upwards. The eigenfrequencies can then be expressed as 

rpj^ 2n + l + 2 



(46) 



with / = 1,2,3... as before. The "radial" magnetic field component can be written as 



= A^%1 + 1) sin'-i X (cosx) ¥^0. 0), (47) 

where A^^ is a normalization constant. The other magnetic field components are, in an 
orthonormal basis. 



Ha 



sin X dx 
imA^^ d r 



sin'+ixC^£+i(cosx) 



sinx dx 



sin'+^xC2S(cosx) 

The time factor exp(— icut) is assumed everywhere. The electric field components Eg and E^, 
not given here, follow from Maxwell's equations. Using these field expressions, we obtain by 



deYim 



Im 



sin^ 



(48) 
(49) 



integrating over the volume the following expression for the total energy: 

2-2'r(2n + 2/ + 3) 



TE 
nl 



vr 



a'\A'''\'l{l + l){2n + l + 2) 



(2n + l)![r(/ + l)]2' 



(50) 



As for the TM modes, we write analogously 



UJ. 



TM 



2n + l + l 



(51) 



with n = 0, 1, 2, ... The electric field components can be written as 



Ea 



^ - A''H(l + 1) sin'-^ xQn '^(cosx) Izm, 
A^*'^ d 



sm^+'xCt'\cosx) 



E, 



sin X dx 

imA^^^ d 



sm'^'xCt'\cosx) 



Imi 



sinx dx 



Y, 



Im 



sin 9 ' 



(52) 
(53) 
(54) 



with corresponding expressions for the transverse components Hg and H^. The total energy 
becomes in this case 

TT _ . 2-2'r(2n + 2/ + 2) 



rpTM 



a'\A'"'\%l + l){2n + l + l) 



(2n)![r(/ + l)P 



(55) 
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• We are now able to discuss how frequency dispersion restricts the modes of oscillations in 
the cavity. As for dispersion relation, we may for a dielectric take a Lorentz (or Sellmeir) 
form as mentioned already in the Introduction, 

where ( is the imaginary frequency and c^o the resonance frequency, the latter being a 
soft frequency cutoff. However, there are no thermodynamic restrictions preventing us from 
assuming that there is a simple sharp cutoff at C = t^O) so let us adopt this simple prescription. 
Moreover, when assuming an ideal metal, we have that e{0) = oo. Our dispersion model 
becomes accordingly 

^«) = { (57) 
[1, C > ^0- 

As for resonance frequency ujo, we shall take the same typical value as mentioned earlier, 

Wo = 1.5 X lO^^rad/s. (58) 

Consider now the TE modes, where the eigenfrequencies are given in dimensional units as 
= (c/a)(2n + / + 2). The lowest mode is obtained for n = 0, / = 1 as {ujQ^)min = Sc/a. 
Let us choose the "radius" of the cavity to be small, 

a = 45 nm (59) 

(this radius is large enough to permit use of macroscopic electromagnetic theory in the 
material). Then, {ujQ^)mm = 2 x 10^^ rad/s. This mode can according to fl58|) not exist in 
the cavity; the permittivity in the walls is simply equal to one. 

In the TM case, we obtain analogously from = (c/a) (2n+Z+l) that (cu™)™*?! = 2c/ a 
for n = 0,1 = 1. This is the lowest possible oscillation mode in the cavity (if we assume 
that there is no restriction coming from dispersion at all). Numerically, (t^(^^'^)mm = 1-33 x 
10^^ rad/s. This oscillation mode can thus exist also in the present case, under the given 
conditions. Our choice of values in Eqs. (158|) and (!59|) has thus managed to make the lowest 
possible mode in the cavity to be the only real one. Our example is somewhat extreme, but 
it serves to demonstrate the important effects of dispersion. 

Equations ( l50l) and ( l55i) permit us to calculate the field energy in each mode. In the 
present case only fl55|) is actual, and it gives for the lowest mode 
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